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' We investigate the importance of unitarity corrections to parton evolution in heavy 

, flavor production at the LHC. The gluon distribution is determined with a fit to 

1 HERA data applying a unified BFKL-DGLAP approach, in which the non-linear 

C D ' evolution is described by the Balitsky-Kovchegov equation. First we estimate bb 

production at CDF and DO. Then, cross sections for heavy quark production 
at various LHC experiments are estimated, tracing the impact of the unitarity 
corrections. 



Oh! 

' 1. Non-linear gluon evolution 



HERA measurements found a steep power-like growth of the gluon den- 
sity with decreasing x which would lead to a violation of unitarity at very 
small X values. The steep growth has to be tamed by gluon rescattering, 
corresponding to non-linear effects in the gluon density evolution. 

The standard framework to determine parton evolution is the collinear 
DGLAP formalism. It works rather well for inclusive quantities but, for 
more exclusive processes, the fc^ -factorization scheme is more appropriate 
because both the longitudinal and transverse components of the gluon mo- 
menta are considered. In this framework, the process-independent quantity 
is the unintegrated gluon distribution, connected to the process-dependent 
hard matrix element via the fcj-factorization theorem. Linear evolution of 
the unintegrated gluon distribution may be described by one of the small x 
evolution equations using the fci-factorization scheme, the BFKL ^ and 
CCFM ^ equations. These equations are based on resummation of large log- 
arithmic pQCD corrections, a" ln'"(l/a;), and are equivalent at the leading 
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logarithmic level. 

The very small x kinematic region is also the regime where the growth 
of the gluon density must be tamed in order to preserve unitarity. Recently, 
a successful description of unitarity corrections to DIS was derived within 
the color dipole formulation of QCD. This is the Balitsky-Kovchcgov (BK) 
equation '^'^ which describes the BFKL evolution of the gluon in a large 
target, including a non-linear term corresponding to gluon recombination 
at high density. 

In our analysis, wc determine the unintegrated gluon distribution 
from the BK equation unified with the DGLAP equation following KMS 
(Kwiecinski, Martin and Stasto) 5, 6, 7,8 ^j^g abbreviation KKMS 

(Kutak, Kwiecinski, Martin and Stasto) for the unified non-linear equa- 
tion. The linear part of this equation is given by the BFKL kernel with 
subleading ln(l/a::) corrections, supplemented by the non-singidar parts of 
the DGLAP splitting functions. Thus rcsummation of both the leading 
InQ^ and ln(l/a:;) terms are achieved. The subleading terms in ln(l/a;) are 
approximated by the so-called consistency constraint and the running cou- 
pling constant. The non-linear part is taken directly from the BK equation, 
ensuring that the unitarity constraints are preserved. One expects that 
this framework provides a more reliable description of the gluon evolution 
at extremely small x, where ln(l/a;) ^ 1 and the unitarity corrections are 
important, than does DGLAP. 

The size of the dense gluon system inside the proton is assumed to be 
R = 2.8 GeV~^, in accord with the diffractive slope, :^ A GeV~^, of 
the elastic J/V' photoproduction cross section at HERA. In this process, 
the impact parameter profile of the proton defines the t dependence of the 
elastic cross section, ~ i?^/2, by Fourier transform. 

2. Constraints from HERA and cross checks at the 
Tevatron 

The initial distribution was obtained by fitting the HERA F2 measure- 
ments using the Monte Carlo CASCADE ^^'^^ for evolution and convo- 
lution with the off-shell matrix elements. The fits were repeated both with 
the standard KMS evolution without the non-linear contribution and with 
extended KMS evolution including the non- linear part. The predicted F2 
is equivalent for both linear and non-linear evolution. 

Next, this constrained gluon density was used to calculate the charm 
structure function F| at HERA and gg bb production at the Tevatron as 
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Figure 1. Bottom production, measured by CDF, is compared to predictions using 
CASCADE with linear and non-linear KKMS evolution, (a) The px distribution of B 
meson decays to J/tj). (b) The azimuthal angle, A<^, distribution of bb pair production 
smeared by the experimental resolution. 



a cross-check of the fit and the evolution formalism. We use rric = 1.4 GeV, 
mb = 4.75 GeV and a renormalization scale in as of = Anig + p^. 
Sudakov Form Factor is included in the calculation. The predicted cross 
section was then compared to both HI Zeus and CDF DO 
measurements respectively. The predictions agree reasonably well with the 
data. 

As an example of these cross-checks, in Fig. 1(a) the cross section for 
B decays to J/tp is shown as a function the J/tp pr The KKMS 

gluon density fits the data both in the linear and non-linear scenarios with 
a comparable accuracy to the NLO collinear approach 

In Fig. 1(b), the azimuthal angle distribution between the b and b 
quarks, Aip, is given. The A(p and bb pr distributions are correlated since 
A(j) < 180° corresponds to higher pair pT- Since the fcr-factorization for- 
mula allows the incoming gluons to have sizable transverse momenta, the 
calculated A<j) distribution agrees very well with the data for A(j) > 60° 
with only smearing due to the experimental resolution. It is interesting to 
note that the inclusion of the Sudakov From Factor improves the descrip- 
tion significantly. For a comparison, we also plotted the result as obtained 
without the Sudakov Form Factor (dotted line). The enhancement of the 
data relative to the calculations at low Acp requires further study. 
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Figure 2. (a) and (b) show bb production as a function of pair pT in the ATLAS/CMS 
(a) and the LHCb acceptance (b). The meson pT distribution in the ALICE accep- 
tance is shown in (c). 

3. Heavy quark production at the LHC 

Since the Tevatron measurements are well described using the unintegrated 
parton densities constrained by HERA and convoluted with the ofF-shell ma- 
trix elements, the same approach may be used for heavy quark production 
at the LHC at e.g. y/s ~ 14 TeV. As discussed previously, heavy quark 
production at this energy is already in the region where saturation effects 
may be relevant. 

We computed heavy quark cross sections for various kinematical regions 
of the LHC. In Fig. 2(a), the bb production cross section is computed within 
the ATLAS and CMS acceptance {pr > 10 GeV and \r]\ < 2.5 for both the 
b and b quarks). In Fig. 2(b), the same cross section is computed within 
the LHCb acceptance where the b quark pT can be measured to 2 GeV for 
1.9 < 77 < 4.9. Similarly, we investigated cc production at ALICE, Fig. 2(c). 
In ALICE, it will be possible to measure the down to px ^ 0.5 GeV in 
\ri\ < 0.9. In the computations the same quark masses and scale was used 
as described in the previous section. 

In all three cases the results of the linear evolution (solid line) and the 
results of the non-linear evolution (dashed line) are very similar. There is no 
significant effect observable for non-linear evolution due to gluon saturation. 
For cc production at ALICE saturation effects have been predicted within 
the GLR approach in the coUinear limit (even with a larger saturation 
radius). This result could not be confirmed with our calculations. 
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The presented results, Fig. 2, suggests that hnear gluon evolution and 
AiT-factorization can safely be applied in the discussed kinematical regions 
of the LHC. 
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